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We show by numerical simulations that a non rotationally symmetric body, whose orientation 
is fixed and whose center of mass can only slide along a rectilinear guide, under the effect of 
inelastic collisions with a surrounding gas of particles, displays directed motion. We present a 
theory which explains how the lack of time reversal induced by the inelasticity of collisions can be 
exploited to generate a steady average drift. In the limit of an heavy ratchet, we derive an effective 
Langevin equation whose parameters depend on the microscopic properties of the system and obtain 
a fairly good quantitative agreement between the theoretical predictions and simulations concerning 
mobility, diffusivity and average velocity. 
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The striking contrast between the simplicity of models 
or experimental setups and the richness and complexity 
of the observed phenomena has contributed to generate 
much interest toward the Physics of Granular Media over 
the past two decades A series of reasons, such as the 
macroscopic nature of grains, their inelastic collisions and 
the lack of true thermodynamic equilibrium represent 
an obstacle to a straightforward application of standard 
methods of statistical mechanics. Dilute granular sys- 
tems, the so-called granular gases today are a priv- 
ileged theoretical and experimental benchmark to test 
the fundaments of kinetic theory and of non-equilibrium 
statistical mechanics in general . 

This paper is inspired to a recent numerical experi- 
ment Q where a Brownian ratchet, i.e. a mechanical de- 
vice able to rectify thermal fluctuations is obtained 
in a non-equilibrium system with energy conserving dy- 
namics. As shown by Van den Broeck, this rectification 
can be obtained by coupling the ratchet to two thermal 
reservoirs at different temperatures without violating the 
Second Principle of Thermodynamics. We underline that 
in order to generate a Brownian ratchet, two symmetries 
must be broken: the time reversal symmetry (detailed 
balance) and rotational invariance of the object. 

We depart from this work by proposing an even simpler 
device, the granular ratchet, which contains the minimal 
ingredients necessary to obtain directed motion. It is 
designed to achieve a non-equilibrium stationary regime 
using the inelasticity and the consequent lack of detailed 
balance @ together with the broken rotational symmetry 
to extract work from a single source. 

The granular ratchet model, sketched in Fig.[l] consists 
of a triangular particle (the ratchet) of mass M , shaped 
as an isosceles triangle with base I and angle opposite to 
the base 26*0 and surrounded by a gas of N disks of diam- 
eter (7 = 1 and mass m = 1. The ratchet can only slide, 
without rotating, along the direction x, perpendicular to 
its base and the whole system is enclosed in a squared 
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FIG. 1: Sketch of the 2D model. The triangle is constrained to 
move only in the x (left/right) direction, while its orientation 
is fixed, i.e. it cannot rotate. Gas particles collide against it 
and occasionally receive energy from an external bath, in the 
form of uncorrelated isotropic random kicks. 



box of side L with periodic boundary conditions. The 
+ 1 particles undergo binary instantaneous collisions 
described by the rule: 



Vi = v', - (1 + ay)cy [(v- - v') • n]n. 



(1) 



where v and v' are the post-coUisional and pre-collisional 
velocities respectively, < 1 is the coefficient of restitu- 
tion for that particular collision, taking value ad if both 
objects are disks or value if the ratchet is involved, n is 
the outward-pointing unit vector normal, in the contact 
point, to the surface of particle z, and Cij is a coefficient 
which takes the value 1/2 if the objects are both disks, 
or the value 1/(1-1- i^n^) if j is the triangle, or the value 
e^/(l -I- e^n^) if i is the triangle, where ~ m/M. Be- 
cause of the constraint the vertical velocity of the ratchet 
is always 0. The collision rule ^ conserves the total 
momentum if i and j are disks, but conserves only the 
x-component of the momentum, when the triangle is in- 
volved. If aij — 1 the total kinetic energy is also con- 
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served. Three possible cases may be considered: (i) a 
pure elastic gas where ad — ar = 1, (ii) a mixed gas 
where = 1 and < 1, (iii) a pure inelastic gas where 
< I and < 1- In both cases (ii) and (iii) an exter- 
nal driving mechanism is needed to attain a stationary 
state and avoid indefinite cooling of the system. Here, 
we use a homogeneous random driving because it is the 
simplest and most studied in theoretical literature 0, S] ■ 
In particular, our simulation implements the non viscous 
version of this thermostat: all disks receive, at a constant 
frequency, independent random Gaussian accelerations 
with zero average and = 1 variance. In experiments 
one may easily reproduce such a thermostat by placing 
the grains upon a horizontal plate vibrating at a high fre- 
quency 0. Since the results are qualitatively the same, 
we will reduce the number of free parameters and restrict 
our simulations to cases (i) and (ii), that is keep = 1 
and vary the inelasticity of the triangle ar only. The sys- 
tem is simulated by means of an event driven Molecular 
Dynamics algorithm. The triangle is "smoothed" by ap- 
proximating its three vertexes with arcs of circle, with 
the condition of tangent continuity along the perimeter. 

The gas is initialized (at t = 0) by assigning to the 
disks non-overlapping random positions and Maxwellian 
velocities with zero average and unitary variance. The 
system forgets its initial configuration and attains a sta- 
tionary state. In the pure elastic case the total energy 
E = TTiivf (where m,; = m if i is a disk and nii = M 
if i is the triangle) is strictly conserved, while in the 
inelastic case, due to the action of the thermostat, it 
reaches a stationary value depending on all the control 
parameters (frequency of the thermostat, collision fre- 
quency, coefficients of restitution, masses, etc.). Numeri- 
cal simulations indicate that the probability distribution 
function (pdf) for the velocity of gas particles and of the 
ratchet are close to a Maxwellian. In the following we 
will indicate as Tg the stationary values of the gas tem- 
perature and Tr the ratchet temperature. All MD results 
given hereafter are obtained using N ~ 1000, L = 500 
(i.e. covered volume fraction ~ 4-10"'^), 6*0 = 7r/6, / = 10, 
giving random uncorrelated kicks to all gas particles ev- 
ery 2000 seconds (gas collision frequency is observed to 
be of the order of 1 collision every 20 seconds). 

We focus here on the statistical behavior of the ratchet, 
whose position and velocity at time t are denoted as X{t) 
and V{t) respectively. Trajectories are averaged over 
1000 realizations starting with different random config- 
urations and discarding the initial transient. Averaged 
trajectories for two particular choices of the parameters 
are displayed in the top frame of Fig. [21 showing that 
when the system is totally elastic no average motion oc- 
curs for the triangle. On the contrary, when inelastic- 
ity is switched on, i.e. ar < 1, even if the external 
driving mechanism acts through random isotropic accel- 
erations and without any privileged direction, the tri- 
angle drifts with average velocity {V) ^ 0. In all our 
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FIG. 2: Top frame: averaged trajectory of the tracer in MD 
(full symbols) and DSMC (empty symbols) with M = 10. 
Bottom frame: rescaled mean squared displacement d2{t) (see 
text) for the same choices of parameters. The power law ~ x 
and ~ are drawn for reader's convenience (dotted straight 
lines) . 



MD simulations we always observed a negative velocity: 
the triangle on average moves toward its base. We also 
studied the tracer self-diffusion, measuring the quantity 
d2{t) = {[d{t)-{d{t))f) with d{t) ^ X{t)~X[0) the dis- 
placement of the tracer with respect to a starting time 
t = 0, taken when the whole system has become sta- 
tionary. This measure is presented in the bottom frame 
of Fig. [21 rescaled following the theory discussed below. 
The usual Brownian behavior with a ballistic first stage 
and diffusive asymptotics is observed. 




FIG. 3: Rescaled average velocities and energies of the tracer, 
both from MD (full symbols) and DSMC (empty symbols) 
simulations. Averages have been obtained with 10'^ MD dy- 
namics and 10* DSMC dynamics. Gas particles have always 
m = 1. In MD, Tg = mlvl) /2 changes when changing the 
parameters (a and M), in DSMC it is always Tg = 1. 



In Fig. [21 the average velocities and energies of the tri- 
angle, measured in MD (gray symbols), are shown. The 
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absolute value of the average velocity of the tracer in- 
creases as elasticity is reduced, showing the origin of 
the effect. As the nice collapse suggests, the ratchet ve- 
locity is proportional to the thermal velocity of the gas 
yj2Tg/'m and to \/M. Therefore, its measure becomes 
very difficult in the limit of large Af, being blurred by 
thermal noise. The bottom frame of Fig. [3] illustrates 
the decrease of the average ratchet energy, with respect 
to Tg, with increasing inelasticity. For smaller ratchet 
masses the effect is stronger. We have also considered 
different values of 9q (not shown in figure): it appears 
that reducing this angle is a way to increase the ratchet 
effect, i.e. increases. This is consistent with the 

fact that non-zero motion of the triangle has origin in its 
asymmetry. Our next step is obtaining some analytical 
predictions to be compared with numerical observation. 

In the dilute gas limit, it is reasonable to study the 
ratchet dynamics by means of a linearized Boltzmann 
equation for its velocity pdf, P(V,t), which can be writ- 
ten as a Master Equation (ME) for a Markov process [1] : 

dP{V, t) 



dt 



dV' [W{V\V')P{V',t) - W{V'\V)P{V,t)] 

(2) 



where the transition rate is: 



W{V\V') 



d9SF{9) I d< / dv'ypcf,{v',,v'y) 

S[V-Vpo,t{V',^.ar,e)]i3) 



(V' -v')-he[{V' ^v') ■ 

with p the density of the gas, Vpost the post-collisional 
ratchet velocity (see eq. (H}), the Heaviside step func- 
tion, S the perimeter length, n = (sinfl, — cos0) and 
for the triangle SF{0) = ^-^{2 sin6'o5(6' - 37r/2) + 
d{9 — 9o) + 6[6 — (tt — 6*0)] }• Following numerical evi- 
dence we approximate the velocity pdf of the gas, 0(v), 
by a Maxwellian with zero mean and variance Tg. It 
is straightforward to verify that detailed balance, in the 
form P{V)W{V'\V) = Pi-V')Wi-V\ - V), holds if 
0(v) is Gaussian and 0^ = 1- 

In order to gain a deeper insight it is convenient to ap- 
proximate the ME by a Fokker-Planck equation (FPE), 
from which we can extract the analytical expression of 
the drift and diffusion terms. This is achieved by express- 
ing the r.h.s. of eq. ([2]) by means of the Kramers-Moyal 
(KM) expansion 



dP{V,t) _ ^ (-1)" 



dt 



d 

dV 



[UV)P{V,t)] (4) 



where jn{V) = J dV'{V' - V)''W{V'\V). By retaining 
only the first two terms we obtain the sought FPE, which 
can be still simplified by expanding these terms in the 
small parameter e. The resulting expressions suggest a 
simple physical picture, which can be illustrated with the 
help of the Langevin equation associated with the FPE: 

V{t) = -^Vit) + ^+Tit) (5) 



The quantities 7, F, T{t) are effective parameters and are 
related to the original parameters by 



7 ~ A-qple 



2ttM 



(1 + sin 6*0) 



F T 
— = -pWil 
M ' M ^ 



sin2 eo)T^{l-rj) 



(6) 



(7) 
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5{t-t') {T(t))=0 (8) 
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(9) 



Hence, for < 1 the ratchet drifts with an average neg- 
ative velocity (V{t)) = F/{AIj). Indeed, the net velocity 
vanishes linearly with e — *■ and is very tiny for massive 
ratchets. It is of interest to observe that in virtue of 
eq. (O the net driving force is proportional to the tem- 
perature difference Tg — Tr, so that the tracer and the 
gas temperatures play role analogous the two reservoir 
temperatures of the Brownian ratchet model. 

To corroborate the robustness of the results we went 
beyond the Maxwellian approximation for the gas pdf 
and used the so-called first Sonine correction, which is 
tantamount to assume 0s(v) = exp ^— (1 -|- 



C2S2imv^ /2Tg)) with S2{x) 



- 2x ■ 



1 hJ. In this 



case one may perform the calculations retrieving a cor- 
rection to formula ([9]) of the form 775 = 77^1^, which 
makes possible the case T^/Tg > 1, i.e. {V) > 0. Note 
that in the literature 77, i.e. the ratio between the gas 
temperature and the tracer temperature, has been calcu- 
lated also for spherical particles in |11|. 

Eqs. (IH])-® indicate that a Fluctuation-Dissipation 
Relation (FDR) holds, in contrast with the small viola- 
tions of FDR reported in studies of different models of 
granular tracers However, the validity of the FDR 

in our case is an effect of the truncation of the KM expan- 
sion and of the small e approximation considered here. 

As anticipated in Figs. [5] and [3] the validity of the ana- 
lytical theory has been also tested against Direct Simula- 
tion Monte Carlo (DSMC [l3|) which enforces the Molec- 
ular Chaos assumption, used to derive Eq. but in 
principle not verified in MD. In addition within DSMC 
it is possible to fix the desired form of (/)(v) at our will, 
while in MD this depends on the control parameters of 
the system. Fig.[3]displays a good agreement between the 
theory and DSMC results for both observables. When 
e ^ 1 the theoretical results deviate from the DSMC re- 
sults. The comparison with MD results is fair, but as 
a,, decreases systematic corrections appear. In particu- 
lar when AI = 1 we observe a nice agreement with DSMC 
results, suggesting that in this case the only source of mis- 
match with the theory is the high value of e rather than 
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the failure of the Molecular Chaos assumption. At higher 
values of e the MD and theoretical results for M(F^) 
match better than the corresponding results for the av- 
erage velocity, because the latter originates from a higher 
order term in the e expansion, and therefore is subject to 
a stronger relative noise. 
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FIG. 4: Rescaled self-correlation functions of the tracer veloc- 
ity with different choices of parameters, elastic and inelastic, 
from MD and DSMC, against rescaled time 'yt. The bold 
dashed line is the theoretical prediction coming from Eq. ([SJ 

The tracer velocity self-correlation C{t) — {{V{t) ~ 
{V)){V{Q) - (y))), displayed in Fig. H together with 
the theoretical prediction Tj./Mexp(— 7t), gives a direct 
measure of 7 and agrees well with the theory for low in- 
elasticity. For larger inelasticities, at large time, C{t) 
shows a fatter tail, which could be the prelude of some 
transport coefficient divergence, very well known for disks 
in 2D [3]. In our opinion this is the explanation for 
the lack of perfect agreement with the estimates of {V). 
Nevertheless, the exponential decay of C{t) is consistent 
with the observed diffusive behavior at large times, as 
shown in Fig. [2l In particular, from Eq. ^ one can pre- 
dict a diffusion coefficient D = Tr/{M-f) and a general 
formula for the mean squared displacement which reads: 
'^^2d'^ = ^ [1 ^ cxp(— z)] where z = jt. This predic- 
tion is fairly confirmed by MD and DSMC simulation, 
see again Fig. [2l bottom frame. 

We are now in a position to draw conclusions and dis- 
cuss perspectives. Our main result is the unveiling of an 
effect unknown, to the best of our knowledge, in granu- 
lar gases: a directed motion driven by undirected fluctu- 
ations, exploiting only the time irreversibility of inelas- 
tic collisions. Experimental verification of such a phe- 
nomenon may easily be achieved: only technical prob- 
lems (such as keeping the tracer in contact only with 
the gas and far from the external driving, in order to 
reduce noise, e.g. constraining it on a suspended guide) 
must be solved. Predictions can be made for different 
ratchets: for instance, we can consider a symmetrically 
shaped tracer, such as a piston or a disk, but with sides 



(say left /right) made of different materials, which corre- 
spond to different inelasticities, ai and 02, respectively. 
The theory predicts (at first order in e) a drift velocity, 
(V) = V27rmkT{u - l){M/m + r] - l)/[4M(it + 1)] with 
M = (1 -I- ai)/(I + 02). Moreover, other kinds of exter- 
nal drivings can be used: a typical setup, for example, 
receives energy from the boundaries. In this case a box 
is vibrated and the tracer should be constrained to be 
in contact with the gas, but free to move on a Id guide. 
Finally, we cannot rule out that similar mechanisms are 
at work also in more packed situations and contribute 
to spontaneous de-mixing phenomena, such as the Brazil 
Nut Problem (H]. 
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